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ABSTRACT 

Wc study the growth of galactic disks in hvc triaxial dark matter (DM) halos. The hales have 
been assembled using the constrained realizations method and evolved from the linear regime using 
cosmological simulations. The 'seed' disks have been inserted at redshift z = 3 and increased in mass 
tenfold over various time periods. 1 ~ 3 Gyr, with the halo responding quasi-adiabatically to this 
process. We follow the dynamical and secular evolution of the disk-halo system and analyze changes 
in the most important parameters, like three-dimensional DM shapes, stellar (disk) and DM (halo) 
radial density profiles, stellar bar development, etc. We find that a growing disk is responsible for 
washing out the halo prolateness and for diluting its flatness over a period of time comparable to the 
disk growth. Moreover, we find that a disk which contributes more to the overall rotation curve in the 
system is also more efficient in axisymmetrizing the halo, without accelerating the halo figure rotation. 
The observational corollary is that the maximal disks probably reside in nearly axisymmetric halos, 
while disks whose rotation is dominated by the halo at all radii arc expected to reside in more prolate 
halos. The halo shape is sensitive to the final disk mass, but is independent of how the seed disk is 
introduced into the system — abruptly or quasi-adiabatically. It is weakly sensitive to the timescale 
of the disk growth. We also expect that the massive disks arc subject to a bar instability, while 
light disks have this instability damped by the halo triaxiality. Implications of these results to the 
cosmological evolution of disks embedded in asymmetric halos are discussed and so are the corollaries 
for the observed fraction of stellar bars. Finally, the halo responds to the stellar bar by developing a 
gravitational wake — a 'ghost' bar of its own which is almost in-phase with that in the disk. 
Subject headings: galaxies: evolution - galaxies: formation - galaxies: halos - galaxies: kinematics 
and dynamics - galaxies: structure - cosmology: dark matter 



1. INTRODUCTION 

A great deal of effort has been invested during the 
last couple of decades in understanding the large-scale 
structure formation in the Universe. While a substan- 
tial progress has been made in following the hierarchical 
buildup of dark matter (DM) halos, the structure de- 
velopment on subgalactic scales, such as the growth of 
disks, bulges, central black holes, etc., is only now be- 
ing investigated. In this paper wc study the buildup of 
stellar disks immersed in evolving DM halos of an arbi- 
trary shape obtained in cosmological simulations and the 
concurrent symbiotic evolution of the disk-halo system. 
Specifically, we analyze the change in the halo triaxial- 
ity to the growing disk and its feedback onto the disk 
evolution. 

Beyond the active and passive adjustments of the halo 
to the growing disk and vice versa, the disk is prone to 
a number of dynamical instabilities, such as the bar in- 
stability, which can alter the evolution of the system in 
the most profound way. A strong and massive stellar bar 
will facilitate the angular momentum and mass redistri- 
bution within the disk and between the disk and the halo 
(e.g., Athanassoula 2003). The efficiency of this process 
is not yet fully understood but is known to depend on 
a number of parameters such as the dispersion velocities 
in the halo, its mass concentration, etc. 

First, from the most general dynamical considerations, 
the coupled dynamical evolution of a disk-halo system is 
expected to depend on the prevailing shape of the DM 



halo and its embedded disk as well as on the ability of 
the disk to develop a strong stellar bar. While DM halos 
form on the average as highly triaxial (i.e., prolate and 
flat) systems in cosmological numerical simulations (e.g.. 
Bullock 2002 for a review), they appear nearly axisym- 
metric (i.e., oblate) in the local Universe (e.g., Rix & 
Zaritsky 1995; Merrificld 2002). What proccss(cs) drive 
this axisymmetrization of the halo? 

Second, the halo radial density profile can have an ef- 
fect on its coupling to the growing disk, especially on the 
efficiency of the angular momentum transfer to the halo. 
It was argued that cuspy halos (i.e., Navarro, Frenk & 
White 1997, hereafter NEW) provide a substantial drag 
on the stellar bars (e.g., Weinberg 1985; but see Sellwood 
2006) and the feedback is capable of leveling off the den- 
sity cusp (Weinberg & Katz 2002; HoUcy-Bockelmann, 
Weinberg & Katz 2005), although this latest point is 
a matter of controversy (Sellwood 2003; McMillan & 
Dchnen 2005). 

Baryons have been shown to modify both the radial 
density profile and the triaxiality of the DM halo — 
dumping baryons within the halo dilutes the triaxial po- 
tential (Dubinski 1994; Kazantzidis et al. 2004) and 
clumpy baryons can level off the density cusp (El-Zant, 
Shlosman & Hoffman 2001). However, the halo triaxial- 
ity (at each radius) can also decrease during the violent 
phases of major mergers as they come from random di- 
rections and during the quiescent phases of slow accretion 
(outside the halo cusp) — this type of halo evolution is 
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actually implemented in the current work. 

The fraction of stellar bars at cosmological distances is 
of course an interesting parameter which can provide us 
with clues to the evolution of galaxies over the Hubble 
time. The original claim that the observed bar fraction 
decreases sharply over rcdshifts z ^ 0.5 (Abraham et al. 
1999) has been questioned first because of the specific 
method of bar detection used in this work (Jogee et al. 
2002a,b). More broadly, analysis of the iTSTGEMSi sur- 
vey of about 1,600 galaxies has shown that the fraction 
of strong bars in the optical stays unchanged to z ~ 1, 
i.e., lookback time of 8 Gyr (Jogee et al. 2004), which 
was confirmed for much smaller samples (Sheth et al. 
2003; Elmegreen et al. 2004). Moreover, Jogee et al. 
have concluded that the distributions of bar lengths and 
their axial ratios is largely preserved as well — with a 
corollary for disk and DM halo properties. Overall, this 
means that since z ~ 1 at least, the galaxy evolution is 
probably not driven by major mergers but entered the 
quiescent phase. 

This result can be coupled to a more general issue 
of the DM halo shape evolution. El-Zant & Shlosman 
(2002) used the method of Liapunov exponents to calcu- 
late the fate of bars in mildly triaxial rigid halos. They 
found that the bars cannot be sustained in these configu- 
rations due to the generated chaos. Bcrcntzcn, Shlosman 
& Jogee (2006) have extended this work to include the 
self-consistent models with live triaxial halos, i.e., flat 
and prolate.^ In all cases, bars which formed as the initial 
disk response to the halo prolateness have dissolved over 
~ Gyr timescale, i.e., few bar rotations. The series of 
models with increasing cuspincss in the halo have shown 
that the chaos generated by bars tumbling in triaxial ha- 
los can also wash out the halo prolateness and thus the 
bar may survive. But it remained unclear whether the 
decrease in the halo triaxiality detected by Berentzen et 
al. is related to the halo relaxation because of the in- 
sertion of a massive axisymmetric disk or it is the result 
of the response in the disk to the halo asymmetry. Our 
present work provides a compelling answer to this ques- 
tion. 

The disk formation within a DM halo has been studied 
before under a long list of simplifying assumptions (e.g., 
Sommer-Larsen et al. 2003; Immeli et al. 2004; Gover- 
nato et al. 2006). Here we focus on the disk- halo inter- 
action and its effects on the developing disk morphology 
and the feedback onto the halo shape. To follow the 
disk- halo evolution in a more realistic environment, we 
monitor the halo assembly in the cosmological framework 
starting with a 2.5(T linear perturbation in the expand- 
ing universe, using the constrained realizations algorithm 
(Hoffman & Ribak 1991) to prescribe the initial condi- 
tions. We insert a small 'seed' stellar disk of about 10% 
of its final mass within the halo which experiences a qui- 
escent phase of its dynamical evolution. The disk then 
grows substantially in size and in mass via self-similar 
mass addition over time periods of ~ 1 — 3 Gyr. We ana- 
lyze the subsequent evolution of such a disk-halo system 
and provide a comparison with simplified versions of this 
system. Those include isolated, axisymmetric and triax- 

Galcixy Evolution from Morphology and SEDs (Rix ct al. 2004) 
^ i.e., b > c and a > b, respectively, where a, b and c are the halo 
principal axes. 



ial halos with non-evolving light and massive disks. 

This paper is structured as following: Section 2 deals 
with the details of numerical modeling. Results are given 
in Section 3 and discussed in Section 4. 

2. NUMERICAL MODELING 

We have used the updated version FTM-4.4 of our hy- 
brid code (e.g., Heller & Shlosman 1994; Heller 1995) 
with the total TV ~ 1.3 x 10*^ particles, of which 1.2M 
represent the collisionlcss DM and lOOK — the coUi- 
sionless baryonic components, respectively. The gravita- 
tional forces have been computed using Dehnen's (2002) 
f alcON force solver, a tree code with mutual cell-cell in- 
teractions and complexity 0(N). It conserves momentum 
exactly and is about 10 times faster than optimally coded 
Barnes & Hut (1986) tree code. The gravitational soft- 
ening of 200 pc to simulate both the evolution of the DM 
component and that of the baryons in the disk has been 
implemented. 

We have introduced the following dimensionless model 
units. The spatial distance unit is taken as r = lOkpc, 
the mass unit is M ~ 10^^ M© and the gravitational 
constant is chosen to be G = 1, which results in a time 
unit of r= (r^/GM)^/^ = 4.7 x lO'^yrs. The dynamical 
timescale, Tdyn within 10 kpc is ~ 4 x 10^ yr for the 
pure DM halo. In these units, the velocity is given in 
208 km s^^. The actual physical units arc mostly used 
here. 

2.1. Initial conditions 

Wc have used the constrained realizations technique 
(Hoffman & Ribak 1991; see Section 2.1.1 below) to con- 
struct a DM halo, and followed its assembly in the ex- 
panding Friedmann universe. At redshift z = 3, we have 
introduced a live 'seed' stellar disk into the halo and grew 
it over the next '^1 — 3 Gyr, depending on the model. 

Vacuum boundary conditions and physical coordinates 
have been used ignoring the cosmological constant term 
and assuming the open CDM (OCDM) model with fig — 
0.3, h = 0.7 and erg = 0.9. Here VIq is the current cos- 
mological density parameter, h is the Hubble constant 
normalized by 100 km s~^ Mpc^^, and erg is the variance 
of the density field convolved with a top-hat window of 
radius 8/i^^ Mpc used to normalize the power spectrum 
(see also Romano-Diaz et al. 2006a, b). Because we are 
interested in the evolution on the sub-galactic scales, the 
difference with the ACDM cosmology is small and our re- 
sults should also stand for the latter. The FTM-4.4 has 
been tested in the cosmological context using the Santa 
Barbara Cluster model (Frenk et al. 1999). Further de- 
tails are given in Romano-Diaz ct al. (2006b). 

2.1.1. The Halo: Constrained Realizations 

The constrained realizations algorithm (Hoffman & 
Ribak 1991) was used by applying linear constraints to 
specify the initial density field in specified locations and 
evaluated with Gaussian smoothing kernels. The width 
of a kernel has been fixed and corresponds to the mass 
scale on which the constraint has been imposed. Within 
the context of the cosmological model and assuming the 
power spectrum of the primordial perturbation field, we 
have generated a constrained realization of the density 
field from a random realization of this field. 
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The simulations start at rcdshift z = 120, correspond- 
ing to T = 0.0127 Gyr, and ends at z = 0.0. The DM 
model has been constructed in such a way as to have a 
set of constraints all located at the origin. The corre- 
sponding nested set of perturbations has been designed 
to collapse off-center. These density constraints consti- 
tute 2.5 — 3.5(7 perturbations, where is the variance 
of the smoothed field, which were imposed on 128'^ cubic 
grid with a side of Ah~^ Mpc. The total mass of the DM 
is ~ 10^^/i~^ Mq and it is embedded in a region corre- 
sponding to a mass of 10^^/i~^ M0 in which the overden- 
sity is zero of unperturbed Friedmann model. The three 
nested off-center perturbations of i5 = 6.0, 4.6 and 3.4, 
corresponding to the mass scale of 9.37 x 10"'^''/i~^ M0, 
1.87 X 10"/i"i Mq and 3.75 x 10"/i-i Mq, are designed 
to collapse by Zcoii = 5.7, 3.8 and 2.2, respectively. Over 
the time of the simulation, the density field has been 
sampled with 500-1,000 time outputs. 

For the reference, the halo mass within 10 kpc is about 
1.5 (or 165K particles), from r = 38, or z w 3 (just prior 
to the seed disk insertion) to the end of the evolution at 
z = (t = 239). This is about two orders of magnitude 
above the halo sampling in cosmological simulations. 

2.1.2. The seed disk 

The initial 'seed' disk is a Miyamoto-Nagai (1975) disk 
described by the gravitational potential: 



where the scalelengths Ad and Bd, and the disk mass 
Mjj are given in Table 1. To include the stellar disk in 
our models, we first determine the main principal axes 
of the triaxial halo, based on its moments of inertia. 
We then introduce the disk with its rotation axis be- 
ing aligned with the halo minor axis. As we show in 
Section 3.1, the halo figure is nearly stagnant over the 
Hubble time. The disk is placed with the center of mass 
position and velocities of the main halo, iteratively mea- 
sured in spheres with radii going from 50 kpc down to 
0.5 kpc. Since the halo is not in equilibrium and the ini- 
tial disk mass is negligible compared to that of the halo, 
we do not require the halo to relax on the introduced 
potential. The disk is inserted at a redshift z = 3, well 
after the last major merger (see Section 3). 

2.1.3. Growing a disk within the live halo 

To simulate the process of disk formation in a simpli- 
fied way, the disk is gradually grown from the seed, by 
adding stellar particles to it, based on the probability 
given by the Miyamoto-Nagai density distribution. The 
central volume density of the disk is kept constant dur- 
ing this time period, and we adjust (i.e., increase) the 
radial and vertical scalelengths of the disk correspond- 
ingly, keeping the ratio Bd/Ad constant. Recent analy- 
sis of the GEMS galaxies supports such a disk buildup, 
when the surface density does (nearly) not change with 
redshift while the disk mass increases with time (Barden 
et al. 2005). 

Since the equatorial plane of the disk may change ow- 
ing to a (slow) inner halo figure tumbling in 3-D, we de- 
termine the instantaneous disk plane from its moments of 



TABLE 1 
Disk model parameters 



Model 


Md 


Au 


Bd 


Rcut 




Dl 


0.0612 


0.132 


0.023 


1.16 


10 000 




0.6120 


0.284 


0.050 


2.50 


100 000 


D2 


0.1838 


0.132 


0.023 


1.16 


10 000 




1.8375 


0.284 


0.050 


2.50 


100 000 



Note. — Columns: (1) the disk model; (2) 
mass; (3) radial scalelength; (4) vertical scale- 
length; (5) outer disk cutoff; (6) number of par- 
ticles. The upper lines in Dl and D2 correspond 
to initial and final values of disk parameters, at ti 
and Tf 



TABLE 2 

Simulations 



RUN Disk Ti [Gyr] Tf [Gyr] Comments 



Note. — Columns: (1) the model; (2) the disk model (see 
Table 1); (3) the start time of the disk growth; (4) the end time 
of the disk growth 



inertia. Then, based on the radial force component at the 
position of the new particles, we assign the correspond- 
ing circular velocity (without any additional dispersion 
or z-velocity component (relative to the disk equatorial 
plane). 

This method, strictly speaking, does not grow the disk 
under exact equilibrium conditions, but we do not expect 
this in a realistic disk formation either, i.e., most likely it 
involves some degree of relaxation as well. We impose a 
linear mass growth over a specified time period, starting 
with rcdshift z = 3 (see Table 2). Overall, the disk mass 
grows by a factor of ten over ^ 1 Gyr or ~ 3 Gyr, which 
corresponds approximately to At = 22 and 60. Because 
the disk center-of-mass position and the 3-D orientation 
of the disk depend on the large-scale motions within the 
surrounding halo and on the streamers on even larger 
scale in the DM, we have developed a new diagnostics to 
address the evolution of disk intrinsic parameters, such 
as the bar pattern speed. 



AO 








pure live DM model 


Al 


Dl 


1.81 


2.84 


default model 


Alb 


Dl 


1.81 


2.84 


as Al, growing frozen disk 


A2 


Dl 


1.81 


2.84 


growing thin disk, live halo 


A3 


Dl 


1.81 


2.84 


growing disk, frozen halo 


A4 


Dl 


1.81 


4.87 


slowly growing disk, live 










halo 


A5 


Dl 






non-growing disk, live halo 


A6 


Dl 






low J disk, live halo 


Bl 


D2 


1.81 


2.84 


more massive disk, live halo 


B5 


D2 






non-growing massive disk, 










live halo 


CI 


Dl 


1.81 


2.84 


growing disk, live axisym- 










metric halo 


C6 


Dl 






non-growing disk, live axi- 



symmetric halo 
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3. RESULTS 

3.1. Evolution of the pure DM halo — model AO 

We evolve the halo model from z — 120 till z = 0, or 
r « 0.013 Gyr and 11.3 Gyr, respectively. We analyze 
the evolution of a diskless halo in the cosmological back- 
ground for further comparison with halos that harbor 
light, massive or growing disk models. 



z z 




a 



Fig. 1. — Evolution of the NFW scale radius Rs, the virial 
radius R^ir, the concentration parameter C = R^ir/Rs and the 
virial-to-total mass ratio M^i^/Mtot, as a function of redshift z 
and cosmological expansion factor a in the pure halo model AO. 
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time [Gyr] 

Fig. 2. — Halo mass in the pure halo model AO within spherical 
radii and as a function of time. The different linos show the mass 
evolution within spheres of constant radii, centered on the main 
halo. The dashed vertical lines mark, from left to right, rcdshifts 
2 = 3, 2 and 1, respectively. 

Initially, the DM expands with the Hubble flow, but 
detaches and starts its collapse by z ^ 15. The epoch 
of major mergers ends around z ~ 7, as given by the 
evolution of the NFW scale radius Rs and the halo virial 
radius'^ Rmr- These parameters experience abrupt in- 
creases during each violent event, as shown in Fig. 1 

^ We define the virial radius iJvir from Mvir = (47r/3)piJ^jj., 
where Mvir is the virialized mass and the current mean density 
p = AcPcriti with Ac being the current virial overdensity and Pcrit 
— the critical density in the universe at this time. 




2 4 6 8 10 



time [Gyr] 

Fig. 3. — Halo shapes based on the isopotential contours as a 
function of time in the model AO. The different colors represent 
the prolateness, t^u = 1 — b/a (upper panel), and flatness, f^u = 
1 — c/a (lower panel) at different radii (black: 5 kpc, red: 10 kpc, 
green: 25 kpc, blue: 50 kpc, and cyan: 100 kpc). 



(e.g., Romano-Diaz et al. 2006a, b). Following the vi- 
olent epoch, Rs remains constant at 7 kpc for the 
rest of the simulations. Rvir grows linearly between the 
mergers and during the quiescent accretion phase which 
follows. The latter phase includes the minor mergers as 
well that occur at all rcdshifts leading to the continuous 
growth of the virialized-to-total halo mass ratio (Fig.^. 
Fig. [2| displays the halo mass growth characterized by 
the mass assembly within the spherical radii of 5 kpc, 
10 kpc, 25 kpc, etc. The overall halo assembly in the 
cosmological background is shown in the Animation Se- 
quence 1. After z ^ 3, the central 50 kpc of the halo 
contains about 650, 000 particles. The virialization pro- 
ceeds from inside out: while the inner 10 kpc virialize 
during the first 1 Gyr, the outer 100 kpc virialize after 
2-3 Gyr. The NFW density profile is established early, 
by z - 10. 

The dimensionless cosmological spin parameter for the 
halo (e.g., Peebles 1969), A = J \E\^^^ /G M^^^, calcu- 
lated within i?vir decreases with time from ~ 0.1 at z 9 
to about 0.013 at z 1, or stays a constant. As ~ 0.1 
within Rs, using notation of Bullock et al. (2001). Here, 
E, J and M are the total energy, angular momentum and 
mass of the halo. These values correspond to the tail of 
the average A in numerical simulations (e.g., Barnes & 
Efstathiou 1987), because we neglect the outer tidal fields 
which become important at lower rcdshifts. 

For our study of a disk evolution within the live halo, 
the halo shape, especially its prolateness, is of a prime 
importance. Fig. |21 shows the halo prolateness e^n = 
1 — b/a and its flatness f^n — l^c/a as a function of 
time, based on a numerically stable ellipse fitting (Halir 
& Flusser 1998) of the isopotential contours, for different 
radii. We note that the isopotentials have a clear advan- 
tage over the isodensity contours, being much less noisy 
(Berentzcn et al. 2006). At the same time the values 
of e$ and /$ reflect the integral properties of the sys- 
tem and have significantly smaller values than those of 
Cp and fp used in the literature and based on fitting the 
isodensity contours. 
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Fig. 4. — Pure DM halo model AO. The figure shows the volume density, in units of 10^ Mq kpc~'^, along the halo principal axes (left 
panel). The right panels show the halo prolateness (upper panel) and its flatness (lower panel). The four frames are taken at different 
redshifts, z = 3, 2, I and 0. The overall halo evolution in the cosmological background is shown in the Animation Sequence 1. 



After about 3 Gyr the halo shape evolves in a quies- 
cent manner, largely in response to the slow accretion. 
The halo prolateness and flatness decrease outwards with 
time because of our neglect of the large-scale tidal field, 
i.e., absence of periodic boundary conditions. The former 
becomes small beyond ~ 25 kpc radius, while the latter 
is significant at all radii. Except for occasional bumps, 
that reflect the halo substructure (i.e., subhalos), both 
parameters change little with time. It is important that 
the inner halo remains prolate and in a quasi-steady state 
in this pure halo model. 

In Fig. 0] we display four snapshots of the halo vol- 
ume density along its principal axes at different redshifts. 
These axes have been determined using moments of iner- 
tia. We also show the radial profiles of halo prolateness 
and flatness. At 2 = 3 (or r = 1.81 Gyr) — the time 
of the seed disk insertion, the halo is prolate, with e$i/ 
just below 0.3 at the center, decreasing to below '^0.1 
outside 100 kpc. By z ^ 2, the halo prolateness is largely 
washed out outside the inner 25 kpc, while halo flatness 
remains significant everywhere. This detail is important 
for understanding the evolution of disks immersed in the 
halo, specifically for the bar instability — the outer halo 
becomes largely axisymmetric in a couple of Gyrs after 
z ~ 3. This happens because of the slow accretion of 
the DM in the absence of the large-scale tidal torques. 
However, it is important that the inner halo, within the 
central 20 kpc, remains significantly prolate and hence 
will interact specifically with the immersed disk. 

3.2. Evolution with growing disks 

3.2.1. Inserting a disk: the halo relaxation effects 

Next, we compare the halo evolution when the disk is 
immersed in a number of different ways. The first pos- 
sibility is when a seed disk is introduced either abruptly 
or more adiabatically at z = 3. The mass of the seed 
disk is small, 10% of its final mass and smaller than the 
halo mass within the disk radius by more than an order 
of magnitude. The disk grows to its final mass within 
~ 1 Gyr, i.e., r ^ 22 or in a less intrusive way over a 
longer period of time, ~ 3 Gyr, i.e., t ~ 60. To verify 
that an abrupt positioning of the disk within the halo 
does not produce any visible effects on the halo evolu- 



tion, even more adiabatic insertion of the seed disk is 
performed as follows. The disk is brought up gradually 
over At = 10 or 20, keeping the disk particle positions 
frozen, but correcting for the disk center-of-mass with re- 
spect to the halo. In these test models, the disk reaches 
its seed mass always at z = 3 (r = 38). At this mo- 
ment we unfreeze the particles and grow the disk tenfold 
thereafter. We perform a number of tests where the basic 
halo parameters are compared between the models (see 
Tables 1 and 2). 

We note that the outer part of the halo itself is not in 
a virial equilibrium, and that the seed disk is dominated 
completely by the halo mass at all radii. The major 
parameters we are concerned with arc the halo flatness, 
and especially its prolateness, e^n- We have run 
a number of test models in order to understand to what 
degree the subsequent evolution of the halo shape de- 
pends upon the method of disk insertion. We started to 
bring up the frozen seed disk at r = 18 ('^ 0.8 Gyr) over 
the time of 0.47 Gyr and at r = 28 (~ 1.3 Gyr) over the 
time of 0.94 Gyr. So by z = 3, the halo feels the full 
seed disk potential, and the disk starts to grow tenfold. 
In all cases the model is evolved for more than 4 Gyr 
thereafter. We find that various tracks, displaying the 
host halo prolateness and flatness in these test models, 
are nearly identical up to z ~ 3, when the disk evolution 
takes over and the curves separate. This means that the 
diverging evolution we observe in some of the models is 
not the result of the halo response to the disk insertion 
but of other factors, such as the disk evolution, discussed 
in the next sections. 

Before we discuss in detail the evolution of the key 
models, we show the halo response to the disk for the 
few most relevant cases: model A5 (a non-growing disk), 
model Al (a standard seed disk growing tenfold), model 
A4 (same as Al but a disk growing over 3 Gyr), model Bl 
(a seed disk growing to a massive disk), and model Alb 
(identical to Al, but the disk is growing and remaining 
in the frozen state) . Fig. [S] displays e$// and f^n taken 
at 5 kpc for these models. The most dramatic decrease 
in the halo prolateness we observe in Bl, where, unlike 
in B5, the disk grows in mass. This model loses all of 
its initial prolateness in the inner halo over the first ~ 
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Fig. 5. — Evolution of the halo prolatcncss (left) and flatness (right) at 5 kpc in different model realizations. The model colors are: 
pure DM halo AO (black), our benchmark model with a growing disk Al (red), slowly growing A4 (green), non-growing A5 (blue), frozen 
growing disk Alb (magenta), massive growing disk Bl (cyan). Vertical dotted lines are 2 = 3 and 2 (left-to-right). 



1 Gyr. The model Al loses most of its prolateness as well, 
while A4 mimics its evolution with a somewhat more 
prolate and flat halo. Interestingly, the model Alb with 
the frozen growing disk has a small but measurably larger 
impact on the halo shape than the live Al disk — to 
be discussed in Section 4. The model A5 with a non- 
growing seed disk leads to a much smaller decrease in 
e^H- The pure DM halo model AO is losing less than 20% 
of its initial prolateness over the first 1.5 Gyr after z = 3, 
i.e., almost three times less than the Al. It experiences 
almost no change in the flatness. Note, not only the mass 
is added in a growing disk but also angular momentum, 
as the new particles possess the maximal allowed J at 
the radii of their insertion. 



B1 




ID 20 10 20 10 20 



R [kpc] R [kpc] R [kpc] 

Fig. 6. — Rotation curves for A5 (time Ti = 1.81 Gyr), Al 
{tj = 2.84 Gyr) and Bl (r/ = 4.87 Gyr) models. Dashed line 
shows the halo contribution, dotted line — the disk contribution, 
and solid line — the total rotation curve. 



The rotation curves for models A5, Al and Bl shown 
in Fig. |B1 correspond to the time when the disk reaches 
its maximal mass. The dynamical importance of the disk 
grows along this sequence — while A5 is halo-dominated 
at all radii, Bl hosts a maximal disk, in the sense of dom- 
inating the mass of the inner region. We note that the 
halo mass is dragged inward in our models because the 
disks are inserted nearly adiabatically and initially have 
a negligible mass compared to that of the halo within the 
same radius. This increase in the halo's central mass con- 
centration is an artifact, because of the lack of clumpy 
baryons in our simulations, and leads to an overall in- 
crease in the halo cuspiness. 

3.3. Evolution with the growing seed disk — model Al 



In the initial response of the seed disk to the triaxial 
halo, a bar-like structure forms immediately and fully 
dissolves ~ 1 Gyr later. At the dissolution time, the 
disk has already reached its full mass. The halo triax- 
iality gets reduced significantly, compared to the pure 
DM model AO, and the halo becomes essentially axisym- 
metric. Because the disk insertion is quasi-adiabatic, the 
DM cusp is not destroyed but becomes even somewhat 
more pronounced, which is purely an artifact of our way 
of introducing the disk (see Section 3.2.1). 

Following the dissolution of the original response to 
the prolate halo, the disk stays slightly oval for a pro- 
longed period of time, in excess of ~ 1 Gyr. For brevity, 
we define the disk as being barred when the amplitude 
of m = 2 mode A2 ^ 0.1, in all models. Using this 
convention, the bar in Al re-appears at ~ 4 Gyr, this 
time as a result of a dynamical bar instability. The de- 
tails of the bar growth and evolution are shown in Fig. [T] 
and the Animation Sequence 2 — the bar experiences 
the vertical buckling instability at r ~ 5.9 Gyr, weak- 
ens and continues to grow secularly, in strength and size, 
for another 2 — 2.5 Gyr thereafter, when it saturates. 
The pattern speed of the m = 2 disturbance in the disk 
is very low initially, but when the bar re-appears, it is 
~ 50 km kpc^^ s^^ . This change can be explained by the 
small initial central mass concentration in the model. On 
the other hand, the fully grown disk adds substantially 
to the central mass and, in addition, drags some of the 
DM inward due to the adiabatic way of insertion. 

The bar in Al reaches its full strength and size only 
after the vertical buckling, z ~ 0.5, but the full disk is 
in place by z ~ 2. What is the reason for this delay 
of ~ 4 Gyr, which is comparable to ^ 10 disk rotations 
within the bar-unstable region of 10 kpc? Typically, the 
strong bar appears after 2-3 disk rotations, e.g., as in the 
axisymmetric models LSl-3 of Berentzen et al. (2006), 
if N is sufficiently large and the discreteness noise is rel- 
atively low. We return to this issue in Section 4. 

The initial halo prolateness triggers a response in the 
disk, with an amplitude, A2, comparable to the halo's 
m — 2 amplitude, i.e., to the degree of its prolateness. 
This response in the disk has a position angle which is 
normal to that of the halo major axis and basically rep- 
resents an oval distortion in the inner disk. Both modes. 
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Fig. 7. — Models with growing disks that develop stellar bars. Upper panels: Evolution of stellar bars (black line) and DM halo (ghost) 
bars (dashed cyan line) m = 2 Fourier amplitudes, A2, integrated within the central 5 kpc. Lower panels: Bar pattern speeds Slj,. The 
vertical dotted lines mark redshifts z = 3 to 2 = 1, respectively. The overall evolution of the disk-halo system for the benchmark model Al 
is shown in the Animation Sequence 2, where projections on the individual planes are displayed. 




Fig. 8. — The DM halo in the model A 1. The figure shows the volume density, in units of 10^ Mq kpc~^, along the halo principal axes 
(left panel). The right panels show the halo prolateness (upper panel) and its flatness (lower panel). The four frames are taken at different 
redshifts, 2 = 3, 2, 1 and 0. 



in the disk and the halo, are stagnant, because the halo's 
spin is negligible, and fade away with time. Also, we 
detect a long-lived spiral structure in the outer disk, in 
this and other prolate models, which was also noted by 
Berentzen ct al. (2006). The existence of this spirals in 
a purely coUisionlcss disk within a prolate halo is in a 
sharp contrast with the disk evolution in axisymmetric 
models. 



As the disk becomes bar unstable, we observe that this 
bar always excites a 'ghost' bar in the DM halo (e.g., 
Athanassoula 2005). This ghost bar is a gravitational 
wake in the halo and not a self-gravitating entity like 
an actual stellar bar. Its shape differs from that of the 
stellar bar - it is much shorter and 'fatter' in all axial 
ratios. The ghost bar would dissolve abruptly if the stel- 
lar bar would disappear. It co-rotates with the stellar 
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bar and has the same pattern speed and nearly the same 
position angle. Aftermath the buckling instability (i.e., 
~ 5.9 Gyr), the stellar bar continues its secular growth, 
and the amplitude of the ghost bar increases as well, al- 
though it always remains well below that of the disk bar. 

Wc have used the Berentzen et al. (2006) method of 
isodensity fitting to estimate the bar size. In Al, it is 
roughly tracking the bar amplitude by dropping from the 
initial 5 kpc (bar semi-major axis) to 3 kpc a 1 Gyr 
later, then increasing monotonically to about 6.5 kpc, be- 
fore the secondary buckling (Martinez- Valpuesta, Shlos- 
man & Heller 2006), and leveling off. 

3.4. Evolution with a growing massive disk — model Bl 

Wc start with a three times more massive seed disk 
than in Al (see Tabic QJ. Wc observed some noticeable 
differences both in the disk and the DM halo evolution as 
a result. The halo prolateness within the central 20 kpc 
drops to near non-existent (see Fig. 5). The initial disk 
response to the halo is very limited in time — the onset 
of the dynamical bar instability cuts it off (Fig.EJ. The 
bar which appears during the disk growth strengthens 
substantially in less than 0.5 Gyr, then decays slightly in 
order to resume its growth for another 1 Gyr. It buckles 
vertically at ~ 3 Gyr, which reduces its strength dramat- 
ically. The bar resumes its secular growth thereafter. For 
rcdshifts z < 0.5, the bar amplitude A2 within the central 
5 kpc shows a slight (but monotonic!) sign of decrease. 

The A2 amplitude (i.e., prolateness) in the halo decays 
initially. The halo ghost bar develops alongside the stel- 
lar bar, experiences a sharp drop in the amplitude during 
the stellar bar buckling, and resumes its secular growth 
subsequently. 

The evolution of the pattern speed of the disk bar dif- 
fers in this model from that of Al as well. The stel- 
lar bar appears tumbling dramatically faster at Qb 
90 km kpc^^ s^^. It decays strongly over the evolution 
time of this model. This behavior characterizes an effi- 
cient J-transfer from the disk to the DM halo over the 
first 6 Gyr of the bar life and a nearly complete absence 
of such an interaction in the next 3 Gyr. 

Furthermore, the bar size follows its A2 variation on 
the average. The bar grows rapidly to above 7 kpc at 
r 2.3 Gyr, then drops to ~ 5 kpc over the next Gyr. 
Interestingly, this drop corresponds to the sudden change 
in the bar growth (see above) and is clearly visible in 
Fig. [71 For the rest of the evolution, the bar size stabi- 
lizes and decays very slowly in the 8-9 kpc range. All 
three major characteristics of the bar, namely its tum- 
bling rate, bar strength and size, correlate — all remain 
nearly constant over a prolonged period of time. 

3.5. Evolution with growing disk in axisymmetric halo 
— model CI 

For a comparison, we have repeated the Al model evo- 
lution in an axisymmetric halo. To create this halo, we 
take a 'snapshot' of AO at z = 3 and redistribute the par- 
ticles azimuthally in the random fashion. We kept the 
original planar radius of particles unchanged and main- 
tained the halo flattening. Afterwards we have truncated 
the halo at 75 kpc. We have added the disk to the model 
and relaxed it for Ar = 50. Finally, we have assigned ve- 
locities in the disk based on the new axisymmetric halo 
particle distribution. 



Naturally, the seed disk-halo system is in dynamical 
equilibrium at the time of insertion, hence no initial disk 
response to the halo is expected (and observed). The bar 
starts to grow in 1 Gyr and buckles subsequently, then 
resumes its growth which saturates secularly (Fig. Cj). 
The ghost bar behaves similarly and weakens during the 
stellar bar buckling, with the follow up growth. 

The stellar bar size grows substantially over the period 
of ~ 2 Gyr. At the buckling time, the bar size stabilizes 
around 6 kpc and then experiences a very shallow secular 
increase to ~ 7.5 kpc, dropping 0.5 kpc at the end of the 
evolution. 

3.6. Evolution with slowly growing disk — model AJ^ 

In this model wc grow the disk over 3 Gyr after its in- 
sertion. The initial response in the disk is the same as in 
Al and so is the final structure of all components (disk, 
bar and the halo). The halo triaxial shape is similar as 
well, with its prolateness remaining slightly larger than in 
Al (Fig. 5). The (initial) halo prolateness is correspond- 
ingly more long-lived. The main difference between A4 
and Al is that the stellar bar starts to form ^ 2 Gyr 
later (Fig. [7)) — which is consistent with the similarly 
slower disk growth in this model. We relate this delay 
in the bar growth to a protracted addition of the angu- 
lar momentum (J) to the disk, as the particles are added 
with the maximal J at each energy. To test this, we have 
run a model with lower J per added particle, model A6 
below. 

3.7. Evolution with a low angular momentum growing 

disk — model A6 

In order to decrease the angular momentum of the 
newly added particles to the growing disk (see model 
A4), we have re-run the Al by assigning only 85% of the 
circular velocity for the new particles. The particles have 
been added only along the major axis of the disk (deter- 
mined from the moment of inertia) in order to provide 
some degree of azimuthal coherency to the growing disk. 
It seems that the disk in A6 is more centrally concen- 
trated than in Al because of the large radial excursions 
of newly added particles with less than maximal J. 

The resulting disk evolution differs from that of Al 
and A4 (Fig. Ej). The stellar bar appears much earlier at 
~ 3 Gyr and its nonlinear growth rate is much higher. As 
a result, the stellar bar buckles substantially earlier than 
in Al. However, asymptotically its evolution converges 
to that of Al, both in strength and in pattern speed. 
The corresponding initial halo response decays faster as 
well and the ghost bar strengthens in tandem with the 
secular growth of the stellar bar. 

Hence, addition of particles on the lower angular mo- 
mentum orbits triggers the bar instability with a higher 
growth rate, which appears to be sensitive to the rela- 
tively mild removal of J, only 15%, from the disk. In 
fact, in A6, the stellar bar starts to form immediately af- 
ter the initial disk response has been washed out (Fig.jTJ 
— about 1 Gyr earlier than in model Al. This confirms 
our previous suggestion that our method of adding mass 
to the disk, i.e., with maximal angular momentum at 
each radius acts as to delay the bar instability. To test 
this, we have run the C6 model, which differs form A6 
only by having an axisymmetric halo. We shall return to 
this issue in Section 4. 
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Fig. 9. — Same as the upper frames of Fig.|7| but for the models with 'failed' stellar bars. 



3.8. Evolution with growing thin disk — model A2 

Instead of giving the newly added particles vertical po- 
sitions based on the Miyamoto-Nagai distribution func- 
tion, we have assigned their positions randomly, between 
z =±100 pc. This was done in order to make the growing 
disk colder and thus more susceptible to the bar instabil- 
ity. Fig. [7| shows that the stellar bar appears somewhat 
earlier, grows faster, but asymptotically is weaker than 
in Al. The halo shape is nearly identical to that of Al. 

3.9. Evolution with growing disk in a frozen triaxial 
halo — model A3 

This model was run in order to test the disk response 
to the undiluted prolatcncss of the rigid halo potential. 
The disk shows a strong initial response which decays 
abruptly to some residual amplitude which is maintained 
by the halo prolateness. The disk grows as in Al but no 
stellar bar appears over the Hubble time of its evolution 
(Fig.O. There can be a dual reason for this: cither the 
mass of the disk is so small that the outer disk cannot 
absorb enough of the angular momentum of the inner 
few kpc, or the actual bar instability is damped by the 
substantially triaxial halo, as argued by El-Zant & Shlos- 
man (2002) and Bcrcntzcn ct al. (2006). Wc have ruled 
out the former possibility by evolving the same disk in a 
frozen axisymmetric halo — a strong bar has appeared in 
this case. Hence this model confirms both the theoretical 
predictions and the previous live numerical simulations 
that the bar instability is strongly damped by the halo 
prolatcncss, especially in frozen halos which lack the abil- 
ity to adjust and to wash out their asymmetry. 

3.10. Evolution with a non-growing disk — model A5 

We start with the default disk, as in the Al model, 
but the seed disk does not grow. As in Al, a bar-like 
feature forms due to the initial response to the triaxial 
halo and dissolves linearly with time. The residual A2 
amplitude both in the disk and in the DM halo have 
the same value (Fig. EJ and can be simply explained as 
driven by the halo prolateness. Because the mass of the 
disk remains small compared to the halo mass within 
the same radius, the disk has no visible effect on the 
halo shape evolution, except that its insertion reduced 
the halo triaxiality somewhat, compared to the pure halo 



model AO. But this effect is substantially smaller than 
for the growing disk in Al. The model A5, in fact, tests 
to what degree the seed disk affects the halo evolution. 

3.11. Evolution with a more massive non-growing disk 
— model B5 

One can suggest, based on the evolution of the Bl, 
that the main reason for the difference in the bar growth 
is the initially more massive seed disk. However this is 
contradicted by the B5 model, where we use the same 
massive seed disk but keep its mass constant. Overall 
impression of B5 is that it closely resembles the A5 model 
(Fig. O. This test model, therefore, provides a strong 
argument that it is not the initial mass of the disk which 
is important but rather its final mass. 

3.12. Bar sizes, disk radial scalelengths and central 

bulges 

In all of the models, bars appear to evolve only in ax- 
isymmetric halos or halos which lost their prolatcncss 
during the disk growth. The timescale of the disk growth 
has no effect on the resulting bar size: models Al and 
A4 develop nearly identical bars with scmimajor axis of 
rf, ^ 8 kpc and a very similar history. Axisymmetric 
disk in CI leads to the same bar as in Al, and so is the 
massive disk in Bl. The low angular momentum disk in 
A6 results in a somewhat smaller bar of r;, ~ 5 — 6 kpc. 
This can be explained by a more centrally concentrated 
disk here. Overall, the bar maximal size appears to be 
close to the DM cusp size Rs ^ 7 kpc in AO. 

The disk growth in our models leads to the increase in 
the central density. Wc find, however, that models that 
host stellar bars develop much larger central densities af- 
ter the period of the bar instability. More massive disks 
growing adiabatically drag in the DM which results in 
substantially higher central densities, especially in Bl. 
Therefore, unlike in Bcrcntzcn ct al. (2006), wc find that 
the surface density in the current disks can be approxi- 
mated by a single-exponential law plus a central 'bulge,' 
instead of a double-exponential law. This 'bulge' region 
corresponds to the axisymmetrized stellar bar part of the 
disk, and wc refer here to the bulge only in the context 
of a density fitting. 

The surface density in the outer disk, beyond the 
bar, was fit with an exponential profile, S(r) ~ 
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I]oCxp(— r/j'd), where So is the normahzation constant 
and rd is the disk exponential scalelength. The result- 
ing exponential law has been extrapolated to the center 
and the excess surface density was considered as a bulge 
contribution. The bulge surface density then has been 
fit by the profile S(r) = Eiexp[— (r/r;,g)^/"], where n — 
the shape parameter corresponds to the pure exponential 
{n = 1) or de Vaucouleurs (n = 4) profiles, ri,g is bulge's 
radial scalelength and its central density (Scrsic 1968; 
Binney & Merrifield 1998). For Al, we find that n « 0.8 
during the entire simulation, and r^g ~ 2.3 kpc when the 
full disk is in place. The scalelength ri,g increases some- 
what to 2.7 kpc until the bar buckles at 5.9 Gyr and then 
decreases secularly towards 2 kpc, due to the action of 
the bar which becomes more centrally concentrated after 
the buckling instability. The ratio r{,/rbg changes, there- 
fore, from ^ 1.9 at the buckling to ~ 4 at the end of 
the simulation. The exponential scalelength of the outer 
disk stays around 3.3 kpc after the disk stops growing. 

3.13. Disk tumbling 

All disk models in our simulations have been permitted 
to move freely in response to the streamers in the assem- 
bling halo and interactions with halo substructure. We 
have detected that the disk plane in model Al starts to 
tilt very slowly around an axis in its plane about 2 Gyr 
after the disk insertion (e.g.. Animation Sequence 2) and 
so in other models as well. This motion has a counter- 
part in the inner halo, within 50 kpc and the equatorial 
planes of the disk and the halo remain aligned over all 
the length of the simulation. The reason for this mo- 
tion in the halo lies in the streamers excited by the halo 
substructure. As a result, we typically do not detect a 
warping of the disk midplane, except for some short peri- 
ods of time. The disk also does not precess for the same 
reason. 

3.14. Velocity anisotropy in the DM halo 

Hansen & Stadel (2006) have claimed a correlation be- 
tween the density slope of the DM halo, a = d\np/d\nr, 
and the anisotropy parameter, f3 = 1 — a'^/a'^, where 
at and cr,. are tangential and radial velocity disper- 
sions in the DM. Specifically, they propose a correlation 
/? = — 0.2(a + 0.8) over about three orders of magnitude 
in r within the DM cusp, with an intrinsic scatter in 
P of 0.05, for slopes with a > —3. We have measured 
both a and /3 for two of our models, pure DM model 
AO and our benchmark model Al. Our best fit for AO 
halo is (3 = — 0.17(q: — 1.68) for a comparable range in 
r (Fig. I10|l . The halo in Al fits the same relation as 
in AO. Hence, we confirm the slope in the correlation 
proposed by Hansen & Stadel, but find the curve being 
significantly offset from their value. 

4. DISCUSSION 

We have investigated the evolution of growing stellar 
disks within assembling asymmetric DM halos, by means 
of high-resolution A^-body simulations. The halos have 
been constructed using constrained realizations in the 
linear regime and evolved from cosmological initial con- 
ditions. We have inserted the 'seed' disks at the center 
of mass of the innermost halo, at the redshift z = 3, after 
the end of the violent (merger) phase in the halo evolu- 
tion. The seed disks, with 10% of their final mass, have 
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Fig. 10. — Inner dispersion velocity anisotropy parameter vs. 
density slope parameter a = dlnp/dlnr for the DM halos in models 
AO (lines) and Al (dots) at different times. The color coding is as 
follows: black at 1.81 Gyr, red at 2.84 Gyr, green at 7.09 Gyr and 
blue at 11.3 Gyr. 



been grown over a period of 1 — 3 Gyr thereafter. These 
disks have been allowed to move freely within the halo 
in response to the interactions with the inner and outer 
substructure and with the streamers within the halo. For 
comparison, we have run additional models with more 
massive disks, non-growing disks, frozen disks, or disks 
in frozen halos, etc. 

Our main result is that a growing disk is responsible for 
washing out the halo prolateness and, to a lesser degree, 
the halo flatness, over a period of time comparable to 
its growth. We find that massive disks which contribute 
more to the overall rotation curve in the system are more 
efficient in removing the prolateness in the halo. The halo 
figure rotation remains negligible during this process. 

In a number of models displayed in Fig. we show 
that the halo shape is very sensitive to the final disk 
mass, while being reasonably independent of how the seed 
disk is introduced into the system — abruptly or quasi- 
adiabatically. The timescale of the subsequent growth 
of the disk mass has a small effect on the halo shape as 
well — such a halo remains slightly more triaxial (model 
A4). On the other hand, when the growing disk was 
kept frozen, the halo lost slightly more of its prolateness 
(model Alb). This has happened because the disk was 
not able to adjust its shape to that of the surrounding 
halo in the latter case. 

The corollary is that the maximal disks most probably 
reside in nearly axisymmetric DM halos. On the other 
hand, light disks whose rotation remains halo-dominated 
at all radii have progressively smaller effect on the halo 
shapes and are expected to reside in somewhat diluted 
but prolate halos. Consequently, we envisage that mas- 
sive disks are prone to development of large-scale stel- 
lar bars (unless additional adverse factors prevail) , while 
light disks have the bar instability damped by the halo 
triaxiality. 

These corollaries from numerical simulations are com- 
plementary and bear interesting observational and evolu- 
tionary consequences for disk galaxies embedded in DM 
halos. Previous theoretical study (El-Zant & Shlosman 
2002) and numerical simulations (Berentzen et al. 2006) 
have argued that stellar bars and triaxial halos are in- 
compatible, by-passing the issue of how such systems can 
form in the first place. Dubinski (1994) has shown that 
growing an analytical spherical potential inside a triaxial 
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halos will reduce the halo asymmetry but residual triax- 
iality will remain. Kazantzidis et al. (2004) have used 
cosmological simulations of the DM and baryonic mat- 
ter to show that cooling within the halos leads to more 
spherical halos than in adiabatic simulations. Models 
presented here resolve the issue of bar-asymmetric halo 
incompatibility by showing explicitly that the growing 
disk is primarily responsible for erasing the halo prolate- 
ness and diluting its flatness, and that long-lived bars 
emerge only in models which become nearly axisymmct- 
ric. 

Substantial differences in the subsequent evolution of 
disk models can be detected, depending on their growth 
history. This is clearly reflected in the growth rate of 
the bar instability. The bar formation is considerably 
delayed in Al and proceeds slower in CI compared to 
lower angular momentum disks in A6 (Fig. 0) ^^nd C6. 
The former disks have grown by a particle addition to the 
circular orbits, i.e., with the maximal angular momentum 
J for the particle orbital energy. As the bar growth is 
associated with the loss of J from the unstable region, 
the addition of J should in principle stabilize the disk, 
temporarily. This, however, does not explain why is it 
that the bar instability is postponed by a number of disk 
rotations after we stop adding mass to the disk in Al 
but appears shortly thereafter in A6. If this difference 
arises solely from the additional 15% of J in Al, one 
should observe the same delay in model CI. However, a 
substantially smaller delay appears in the (axisymmetric 
halo) CI. 

We note that the original halo prolateness is preserved 
for a longer period of time in Al than in A6, as the direct 
comparison of halo's m = 2 amplitudes in Fig. shows. 
One observes in this Figure that the stellar bar starts 
to grow when the halo's A2 amplitude drops below 0.1, 
i.e., small e^// — in tandem with calculations based on 
the Liapunov exponents and in self-consistent A^-body 
simulations of disk evolution in live triaxial halos (El- 
Zant & Shlosman 2002; Berentzen et al. 2006). 

We have run a test model (C6) with an axisymmet- 
ric halo, as CI, but where particles have been added 
with 85% of J for the circular orbits — similarly to A6. 
Its stellar bar appears immediately after the disk has 
reached it final mass, by z 2. What is evident is that, 
in both A6 and C6 models, the stellar bar pattern speed 
has a maximum which is clearly larger than in Al and 
CI. The simple explanation of this is that models with 
lower J, A6 and C6, are more centrally concentrated than 
the latter models and their central dynamical timescale 
is accordingly shorter. Hence, at least in the linear stage, 
the growth rate of the dynamical instability is larger in 
A6 and C6. It appears that a combination of an addi- 
tion of J and different central concentration govern the 
bar/disk evolution in the models described above. 

All models developing stellar bars also develop ghost 
bars in the halo that are nearly aligned with the disk 
bars. Interaction between these bars is a major contrib- 
utor to the slowdown of the stellar bar. The ghost bar 
acts as a gravitational wake and its orbital structure is 
fundamentally different from that of the stellar bar — in 
particular, the DM particles arc not trapped within the 
ghost bar as they are within the stellar bar. Ghost bars 
are sensitive enough to respond to the vertical buckling 
in the stellar bars and they grow in tandem with the 



stellar bars, as seen in Fig. |7| However, they should 
not be confused with the original halo prolateness which 
also contributes to the m = 2 amplitude, A2, in the halo 
for the first few Gyrs. We find that this latter contribu- 
tion serves as a reliable indicator for the growth/decay of 
the stellar bars — those grow when halo's A2 amplitude 
drops below 0.1. 

The density profile in the disk can be fit by a 
single-exponential profile plus a central bulge which is 
somewhat steeper than an exponential. This differs 
from double-exponential disks that have developed in 
Berentzen et al. (2006) models. The reason for this 
difference appears to be the quasi-adiabatic insertion of 
the disk in the current work. This quiescent way of in- 
sertion drags in the DM, without destroying the central 
cusp. Berentzen et al. (1998) has found a similar bulge 
shape for pure coUisionless models, which has been mod- 
ified with the addition of the gas component in the disk. 
Overall we find that disk-to-bulge scalelength ratio stays 
in the range of 1.2 — 1.6, while bar size-to-bulge scale- 
length ratio increases from about 1.9 to 4. In all models, 
the bar has been roughly confined to within the NFW 
scalelength Rs- 

We find that the DM halo cusps are preserved in our 
numerical simulations, although the quality of the NFW 
fit is significantly degraded with the growth of the disks. 
The reason for this apparent robustness of the cusps lies 
in the quasi-adiabatic insertion of our seed disks which, 
in fact, reinforce the cusps by dragging in the DM. Even 
in models with strong stellar bars, the cusps survive. It 
is possible that in reality they are destroyed by the ac- 
creting clumpy baryons (e.g., El-Zant et al. 2001, 2004) 
or by other similar processes (e.g., Tonini, Lapi & Salucci 
2006) prior to the disk formation. This issue is beyond 
the scope of this work. We also did not address the 
possible effect that the halo cuspiness has on the halo's 
triaxiality (e.g., Berentzen et al. 2006). 

How does the evolution of disk-halo systems described 
in this work fit within the broader issue of disk galaxy 
evolution within the cosmological context? While numer- 
ical simulations have shown that halos form as prolate 
and fiattened entities (e.g.. Bullock 2002), they appear as 
nearly axisymmetric (oblate) in the local universe (e.g., 
Merrifield 2002). On the other hand, the high frequency 
of stellar bars in the local and intermediate redshift uni- 
verse (Jogee et al. 2004; Elmegreen et al. 2004; Sheth et 
al. 2003) poses a potential problem because bars and pro- 
late halos are incompatible (El-Zant & Shlosman 2002; 
Berentzen et al. 2006). A process which washes out the 
halo prolateness, therefore, can be, in principle, responsi- 
ble for the present high fraction of barred disks. For this 
to happen, the reduction of a halo triaxiality should be 
achieved quasi- adiahatically. Otherwise, a galactic disk 
will be heated up and stabilized against the bar forma- 
tion for a prolonged period of time. Major mergers and 
even frequent minor mergers can lead to such an 'over- 
heating' of stellar disks. Furthermore, the major mergers 
will typically destroy the disks. 

A gradual buildup of a galactic disk within an ini- 
tially asymmetric DM halo can play the role of the re- 
quired quasi-adiabatic process which resolves both issues 
at once, i.e., the halo prolateness and the high observed 
bar fraction. Jogee et al. (2004) have argued in favor 
of such a quiescent disk evolution over the last 8 Gyrs, 
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based on the inferred optical fraction of bars and on the 
distribution of bar ellipticities and sizes over z ^ 0.2 — 1 
in recent GEMS survey with the RST. Moreover, Barden 
et al. (2005) have shown that the disk mass grows faster 
with redshift than its surface density. This requires a 
gradual increase in the disk radial scalelength with z — 
a method implemented here to grow the seed disks. 

The scenario that emerges from this gradual buildup of 
disks within the triaxial DM halos, in principle, can ex- 
plain the current population of stellar bars in disk galax- 
ies. Within this hierarchical merging framework, the DM 
halos of arbitrary shapes form first during the major 
merger epoch. If disks form during this early stage, they 
are probably destroyed shortly thereafter. Disks that 
grow during the next quiescent epoch can be responsi- 
ble for the washing out of the halo prolatcness and, to a 
lesser degree, its flatness. The efficiency of this process 
depends on the increasing disk-to-halo mass ratio. Large- 
scale stellar bars will form under these conditions with a 
delay — first bars would form in systems with more mas- 
sive disks and, therefore, more oblate halos. For lower 
redshifts, progressively more disk galaxies will be sus- 
ceptible to an intrinsic (spontaneous) bar instability as 
contrasted by the tidally-induced bars. 

Hence, bars which are observed in the local universe 
can still be the first generation which formed sponta- 
neously after the host halo prolateness has been washed 
out. Within this scenario, the intrinsic bar formation is 
spread out in time. Thus, the observed fraction of stellar 
bars, and distributions of their strength and sizes reflects 
a steady state between the newly forming bars in more 
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axisymmctric halos and bars destroyed by various pro- 
cesses. The details of the above formation and destruc- 
tion processes are still sufficiently unclear. For example, 
we do not know what is the characteristic timescale of a 
disk buildup. This work explores the range of 1 — 3 Gyrs, 
but, in principle, it can be much longer, ~ 8 — 10 Gyrs, 
as hinted by the results from GEMS, and, therefore, can 
extend beyond z ~ 2. 

To summarize, wc have developed a method to insert 
seed disks and to grow them subsequently in cosniologi- 
cally evolving halos in a quasi-adiabatical way. We have 
investigated the response of the halo shape to the growing 
stellar disk and the feedback of the halo onto the subse- 
quent evolution of this disk. We find that the growing 
disk dilutes the halo triaxiality depending on the final 
disk mass. Maximal disks, therefore, are expected to re- 
side in oblate halos, while disks dominated by the halo 
at all radii can reside in prolate halos. Consequently, we 
expect that the massive disks are prone to the bar insta- 
bility, while light disks have this instability damped by 
the halo triaxiality. 
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